18.785: Analytic Number Theory, MIT, spring 2007 (K.S. Kedlaya) 

The Selberg sieve 



1 Review of notation 

Let / : N ^ C be an arithmetic function, and suppose we want to estimate the sum of / 
over primes. More precisely, let P be a set of primes, and put 

p{z)= n p- 

P<z,peP 

If we define 

S(x,z)= Yl /(^)' 

n<x,{n,P(z))=l 
n<x,n=0 (mod d) 

(with the dependence on P and / suppressed from the notation), we have 

S{x,z) = E fi{d)Ad{x). 

d\P{z) 

Let g{d) be a multiplicative function with 

5(p)G[0,l) (peP); 
g{p)=Q ip^P), 



and write 
Then 



^d{x) = g{d)x + rd{x). 

S{x, z) = V{z)x + R{x, z) 
V{z)^ Uil -g{p)) 

p\P{z) 

R{x,z) = rd{x). 

d\P{z) 

2 The Selberg upper bound sieve 

In the previous unit, we used the combinatorial sieve to construct an arithmetic function 
A+ : N ^ M such that 

A+(l) = 1 
EA+(d)>0 (n>l). 

d\n 



1 



By setting 

d\P{z) 

R+{x,z)= J2 A+(c?)rd(x), 

d\P{z) 

we were able to obtain the bound 

V-(z)x + R-{x,z) < S{x,z) < V+{z)x + R+{x,z), (1) 

but controlling and was rather painful. 

Selberg proposed instead to construct an arithmetic function p : N — > M with = 1 
and 

d\n \ d\n 

In other words, let p be any arithmetic function with p(l) = 1, and put 

A+(n) = Yl Pidi)p{d2)- 

di ,(^2 :lcm((ii ,d2 )=n 

We will typically want \'^{d) = for d> y, for some prespecified number y; to enforce this, 
we may insist that p(n) = for n > ^/y. We call the resulting A"*" an L^-sieve of level y, or 
more commonly a Selberg (upper bound) sieve of level y. 

Let us drop x from consideration by agreeing to only consider functions / with finite 
support. (That is, we replace / by the function vanishing above x.) If we again set 

S{z)^ Y fin) 

{n,P{z))=l 

V^{z) = Y ^^(^)^(^) 

d\P{z) 

= Y p{di)p{d2)g(\.cm{di,d2)) 

di,d2\P(^) 

R-^{z) = >^^(d)ra{x) 

d\Piz) 

= Y P((^l)Pid2)ncm{dud2){x), 
di,d2\P{z) 

we again have 

S{z) <V+{z)x + R+{z). (2) 
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Ignoring the error term R'^{z) for the moment, one can ask about optimizing the main term 
V'^{z)x in the bound (2). This amounts to viewing V'^{z) as a quadratic form and then 
minimizing it. 

For simphcity, we will assume that g{p) G (0, 1) for p E P, and g{p) = for p ^ P. 
(Before we only wanted g{p) e [0, 1) for p G P, but there is no harm in adding to P those 
primes p for which g{p) — into P.) Let h he a. multiplicative function with 

AW- 



1 - g{p) ' 

We can then diagonalize the quadratic form as follows: first, put c = gcd(di, ^2), ^ — di/c, 
b — d2/c to obtain 

^^(^) = E p{ac)p{bc)g{abc) 

a,b,c:abc\P{z) 

- E 9(0)-' E {9{ac)p{ac)){g{bc)pibc)). 

c\P{z) a,b:abc\P{z) 

Note that since P{z) is squarefree, the condition abc\P{z) forces gcd(a, 6) = 1. We now 
perform inclusion-exclusion on gcd(a, b) to obtain 

^^(^) = E E E g{cdm)p{cdm) 

c\P(z) d\P(z)/c \m\P{z)/{cd) 
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m) 



= E 3{c)-' E /^(^) E 3{m)p{ 

c\P{z) d\P{z)/c \m\P{z):m=Q{cd) 

We next substitute e, //e in for c, 0?, and reorder the sum: 



^^(^) = E E/^(//^)^(^)"' E 9{m)p{m) 

f\P{z) e\f \m\P{z):m=0{f) 

E M/)-' ( E 9{m)p{m) 

/|P(z) \m\Pizy.m=Oif) 



Let's put 



so that we have 



^{d) = p{d) g{m)p{m), 

m\P{zy.m=0 (d) 

V^z) = Yl h{dr%df. 

d\P{z) 



Before we can minimize this quadratic form, we must first reexpress in terms of ^ the con- 
ditions we imposed on p. Namely, by Mobius inversion, 

' d\P{z):d=Q{n) 

SO the condition p(l) = 1 is equivalent to 

d\P(z) 

and the condition p{d) — Q lor d > is equivalent to 

e(d) = o {d>^). 

That is, is restricted to a hyperplane. 

Here's where the part comes in. By the Cauchy-Schwartz inequality, 

V+{z) > H-\ E h{d) 

d<^,d\P{z) 

and equality holds for 

^{d)^h{d)H-' {d<^). 

Backing up, we get 

p{d) = p{d)^H-' Yl Mn)- 

^ n<^/d:gcd{d,n)=l 

Putting this together, we obtain the following. 

Theorem 1 (Selberg). Let / : N ^ ]R>o be an arithmetic function with finite support. Let 
P be a set of primes, and put P{z) = Y[p<zp&pP- -^'^^ d\P{^)! write 

Ad^ f{n)^g{d)X + ra{z) 

for X > Q and g a multiplicative function with < g{p) < 1 for all p & P. Let h{d) he a 
multiplicative function with h{p) — g{p){^ — g{p))~^ for all p & P, and put 

H= Yl ^(^) 

d<^,d\P{z) 



for some y > 1. Then 



S{z) = Y ^ ^H-' + E >^''id)r,{z), (3) 

(n,P(z))=l d\P{z) 



for 

A+(n) = J2 P(di)p{d2) 

di ,d2 ;lcm(di ,d2)=n 

p{d) ^ „{d)^H-' M^)- 

' n<VI^/<i:gcd(d,n)=l 

As a somewhat miraculous corollary (due to van Lint and Richert), we obtain 

< n{d)p{d) < 1 (4) 
(exercise); this makes it easy to estimate the error term in (3), e.g., by 

|A+(d)| < (i(^''«5 3)/(log2) 

(exercise) . 

Exercises 

1. Prove (4). (Hint: group terms in the definition of H according to the common divisor 
of d with some fixed number e.) 

2. Deduce (5) from (4), by proving that |A"'"(d)| < 3'^^''^ for u{d) equal to the number of 
prime factors of d. 

3. In the Selberg sieve, prove that if we extend to a completely multiplicative function, 
then 

H>J2 9{n). 



4. Prove that for some c > 0, 



— >ciog^x {x>i). 



n 

n<x 



(Hint: an elementary proof is possible, but one can also use analytic arguments on the 
Dirichlet series C^(s)/C(2s) = En=i 2^^^"^ W-) 

5. Let d{n) denote the number of divisors of the positive integer n. Prove that 

d{n) ~ X log X. 

n<x 

(This is needed for the next problem.) 



6. Use the Selberg sieve to prove that the number of twin primes p < x is 0{x/ log^ x). 
(Hint: put f{n) = 1 if n = m(m + 2) for some m and f{n) — otherwise, then apply 
the Selberg sieve with z — x^^^. You may need some of the earlier exercises as well.) 

7. (Brun-Titchmarsh theorem) Prove that for any e > 0, there exists xq — xo{e) with 
the following property: for any positive integers m, N with gcd(m, N) — 1, and any 
X > max{A^, xo{e)}, the number of primes p < x with p = m (mod N) is at most 

(2 + e)x 



(f){N)\og{2x/N)' 

This is one of several problems in which the Selberg sieve applies to give you a result 
which is off by a factor of 2 from the expected best result. 

8. Prove that 

n<x ^ ' 

then deduce by partial summation that 

(Hint: first prove that the sum X]„ 1/(?T'7(?T')) converges, where ^{n) — W^^^'P-) 

9. Use the previous two exercises to deduce that 

Y,d{p-\)^0{x\ 
where d{n) denotes the number of divisors of n. 



6 



